THE CAPILLARITY PROBLEM FOR COMPRESSIBLE LIQUIDS 
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ABSTRACT. In this paper we study existence and regularity of solutions to the capillarity problem for 
compressible liquids in a tube. We introduce an appropriate space of functions of bounded variation, 
in which the energy functional recently introduced by Robert Finn can be defined. We prove existence 
of a locally Lipschitz minimiser in this class. 



1. Introduction 

Extensive work has been published on the behaviour of capillary (liquid-air or liquid-liquid) 
interfaces when the liquid is assumed to be incompressible. As an authoritative introduction we 
refer to [6] by Finn. Two major approaches have been used to obtain existence and regularity results: 
classical PDE techniques for surfaces of prescribed mean curvature (see, for example, [11, 15, 16]), 
and the functions of bounded variation and sets of finite perimeter setting for minimising the energy 
(see, for example, [4, 10, 12, 14, 15]). 

Results concerned with compressible liquids are very recent and comparatively few, the model 
having been introduced by Finn in 2001 [7], see also [8]. Following that paper we consider a 
capillary tube of cross section Q, C W 1 , which is simply connected and has Lipschitz boundary 
S := dQ. We also assume that it satisfies an interior sphere condition of radius R. 

The capillary surface S is given as a graph of a function u over the domain S7. We assume uniform 
downwards gravity g and consider a compressible fluid of density p. (In the incompressible fluid 
case, p is constant.) 

One can assume prescribed mass M, but the results in the present paper are for an infinite con- 
tainer. 

We consider the energy for a capillary surface to consist of the following components: 
Energy of the free surface (surface tension): 



E s = — I <S>(u;p ) V / l + \Du\ 2 dx; 
Po Ja 



Potential energy: 
Wetting energy: 



W = g / h$(h;po) dhdx; 
Jn Jo 

Ex = —a (3 I $(h;p ) dhds; 
Jt, Jo 



here f3 6 is the relative adhesion coefficient, satisfying \f3\ < I — a with a > 0; a and g 

are the surface tension and gravitational constants; $(/i;po) is the density function depending on 
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height h and pressure p, which we assume to be given by one of the two models proposed by Finn 
[1, 7]. In the following, po and po will denote pressure and density at a reference level u = 0. 

Mass: In the case of a mass constraint, a term AM is added to the energy, where A is a Lagrange 
multiplier and the mass is 

M= / <&(h;p )dhdx. 
Jn Jo 

The total energy (and in particular the wetting energy) need not be positive. 
A smooth minimizer of the total energy Es + W + Ey, + AM will satisfy the Euler-Lagrange 
equation 

(1.1) div ; = u-\ — ——■ — ; + A — on U, 

tJI + \Du\ 2 o- <S>(u;p ) y/l + \Du\ 2 a 

with boundary condition 

1 Du-v 
P = ; on S, 

PO y/l + \Du\* 

via standard calculus of variations techniques. 

The present paper is based on one of the models proposed by Finn for an isothermal fluid: the 
density is assumed to be linear in the pressure, from which one obtains that $(h;po) = poe~ X9h , 
for some positive constant x- 

We may assume that x = h 9 = 1> Po = 1> & = 1; other values of these constants correspond to 
different weightings on the components of the energy (that is, our energy becomes 71-E5 + +72VF+ 
73E-E + 74AM for 7i > 0), and a scaling of the domain Q. 

Then 

Q(h;po) = e~\ 

and the diverse components of the energy are as follows: 

E s = I e-Vl + |ZM 2 dx, 
Jn 

W= I I he~ h dhdx= [ [1 -e- u {l + u)] dx, 
Jn Jo Jn 

E^= [ (3 [ e' h dhds = - [ (3 (l - e~ u ) ds, 
Jt. Jo Jt. 

M= [ [ e~ h dhdx= I (l-e- u )dx. 
Jn Jo Jn 

As we are dealing with the case of an infinite container, we choose A = 0. Without loss of generality 
(this will be shown when necessary, in Lemma 3.6) we may set 7$ = 1, and then seek to minimize 
the energy 

J(u) := E s (u) + W{u) + Ez(u). 

The following results have been recently obtained for the capillarity problem of a compressible 
fluid. 

For slightly compressible fluids Finn [7] introduced the model we are using here. In the case of 
a tube closed at the bottom he found the necessary condition on the mass for existence of a solution 

is M < p \n\/ X g. 
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For a circular tube, Finn and Luli [9] show that for any boundary contact angle 7 with < 7 < ir 
there is at least one symmetric solution of the problem, and that the height of this solution will lie 
above any prescribed level if M is sufficiently large. If 7 < tt/2, the solution is unique among 
symmetric solutions with that mass. 

Finn and Athanassenas [1] follow the classical PDE approach. They include the situation where 

on the right hand side of the prescribed mean curvature equation (1.1) the term ^-u is replaced 

a 

, Po — XPo, X gu _ an( j stU( jy non-constrained case. The results vary depending on the 

regularity of the boundary of the domain: they obtain height and gradient estimates and existence 

of smooth solutions for smooth domains, but only variational solutions for domains with Lipschitz 

boundaries. As with the incompressible case, they observe that existence of solutions depends 

discontinuously on the opening angle of the corners of the domain. In the case with the alternative 

right hand side of (1.1), they show non-existence whenever the domain is small, that is, when 

Po ~ XPo \n\ > -\T,\/3 (here is taken to be constant). 
o-X 

In the present paper we use functions of bounded variation techniques. 

In Section 2, we introduce BV, the space of functions of bounded variation. After a transfor- 
mation of u, the weighted surface area term is well defined in BV. Transforming the remaining 
components of the energy gives us a new energy, J\. 

In Section 3, we prove height estimates. In Lemma 3. 1 and Lemma 3.5 we have two Stampacchia 
type results needed in our case. 

In Section 4, we show that the energy functional is bounded from below, and that a minimising 
sequence for the energy functional is uniformly bounded in the BV-norm. Existence then follows 
via the standard compactness theorem and by the lower semicontinuity of the functional. 

Finally, in Section 5 we show that there exists a locally Lipschitz minimiser. 
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2. The energy in the isothermal case 
As in [13] the BV-seminorms are: 

/ a/ 1 + \Dv^ = sup { / g n+ i + u div„ g dx : g { G C\ (fi) Mi = 1, . . . , n + 1, V] 9i 2 < 1 > 
Jn [Jn ~J J 

and 

/ \Du\= sup I / u div n g dx : gi £ C\ (O) Mi = 1, . . . , n, V] gj 2 < 1 } , 
Jn [Jn ~l J 

where div n is the divergence of the first n components, div n g = J2?=i ^i9i- 

One then defines the spaces BV(fi) := {u G L l (VL) : f n a/I + \Du\ 2 < 00}, and BV + (f2) := 
{u G BV(ri) : u > almost everywhere in Q}. 
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In the case that u G C 1 ^), the surface energy term E$ = f n e~ v + \Du\ 2 dx may be sim- 
plified by writing v = e~ u . Then we can rewrite it as f n \/v 2 + \Dv\ 2 dx, which bears close resem- 
blance to the integral investigated by Bemelmans and Dierkes [2], which was J n \AT+ \Dv\ 2 / 4 dx; 
see also [3]. 

The focus of our investigation now shifts to v, rather than u itself. 
Define 

/ ^v 2 + \Dv\ 2 : = sup < / v (g n+1 + div n g) dx : ft G Cq(Q) Mi = 1, . . . , n + 1, V] ft 2 < 1 \ 
Jn [Jn ~~y J 

Lemma 2.1. Ifv is smooth, J n \fv 2 + \Dv\ 2 = J Q t/^^JDvI 2 dx. 

-| icy 

Proof. We consider the test function g 6 = \e [v 2 + I -CM 2 ] (~Dv, v), where \e is a sequence 
of Cq°(Q) functions with \e < L converging to xn, the characteristic function of ft, in L 1 . Then 

/ ^v 2 + \Dv\ 2 > [ v(g E n+1 +div n g £ )dx 
Jn Jn 

= / g £ n+1 v - Dv g £ dx 
Jn 

= / Xe \/v 2 ~+JDv^ dx 
Jn 

\/v 2 ~T]Dvf 2 dx 

as e — > 0. The other direction is similar. □ 

We note the following fact: 
Lemma 2.2. The quantity i/t> 2 + \Dv\ 2 is finite exactly when v is in BV(fi). 
Proof. Suppose that v is in BV(S1). Then 



/ 

Jn 



/ ^/v 2 T]Dv\ 2 = sup I / v (g n+1 + div n 3) cfe : ft G Cq (U), ft 2 < 1 > 
Jn [Jn ~[ J 

< sup jy (|u| + l)5 n+ i +wdiv„5f(ix : ft G Cq(^), ^Z^ 2 < l| 

< sup 1^ +wdiv„5(ia: : ft G Cq(^), Yl 9 ? - 1 j 

+ sup|y ftj+i|w|dx : 5 n+ i G Co(^), fti+i 2 < l| 



/ 



in 

< 00. 



y 7 ! + \Dv\ 2 + ||u|| L i 



(O) 



On the other hand, if J n \/v 2 + \Dv\ 2 < 00, then v G since if not, we can take ft = for 

i < n + 1 and ftj + i = Xe (where Xe is as m Lemma 2.1) so that 



/ \/v 2 + |Dw| 2 > / |t;|xed2; 



00 
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as e — > 0, contradicting our assumption. Finally we can check 

/ v 7 ! + \ Dv \ 2 = SU P { I 9n+i + v div n g + vg n+1 - vg n+ i dx : g { G Cq (ft), V] s>j 2 < 1 > 
J si [J si ~i J 

< sup | ^ f div n g + v# n+1 : g { G Co(fi), < 1 j 

+ sup j^sf„ + i(l : G Co(ft),fii„ + i 2 < 1 j 



< [ V^ 2 

Jsi 



j 2 + \Dv\ 2 + \\v\\ L i (n) + \n\ 

< oo. 

□ 



Corollary 2.3. If, in addition to the above conditions, v is in BV(fJ), and Vk is a mollification of v, 
then 



[ Jv 2 + \Dv k \ 2 - / V^ + I^l 2 - 
7n v J si 



This may be proved in the same manner as Lemma A 1 of [10]. 

Under the transformation v = e~ u , the wetting energy is = — J* s j3[l — e~ u ]ds = — J* s (3[\ — 
v]ds, where we consider u| s as a trace of v. As in [13, Theorem 2.10] if £ is Lipschitz, each 
function in BV($7) has a trace in L 1 (£). Furthermore, if S also satisfies an interior sphere condition 
with radius i?, then the following estimate holds (see [10], Remark 2): 

(2.1) / \v\ds< / \Dv\ + c R / \v\dx, 

7s 7n 7n 

where cr depends on n, R, and S. 

The integrand of VF, the potential energy term, becomes In /i d/i and so the complete energy, 

in the isothermal case, is 

(2.2) Ji(v)= [ Jv 2 + \Dv\ 2 + f [ Inhdhdx- [ 0[l-v]ds. 

Jsi J si 7i 7s 

Here we are reminded of the energy studied by Claus Gerhardt in [10], which was 

[ Jl + \Dv\ 2 + f [ H(x,h)dhdx - [ f3vds 
J si J si Jo 7s 

dH 

for (3 € L°°(S) and H satisfying the conditions (a) — — > 0, and (b) H(x, h ) > (1 + c), 

oh 

H(x, —ho) < — (1 + c) for some /i > and a given c. The conditions imply that for large values 
of \v\, the potential energy term (the one involving H) is strictly positive and increasing at least 
linearly in \v\. The current case is an improvement on this situation: a strictly positive potential 
energy, which increases like v In v for large values of v. 
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3. Height bounds 

In this section, we assume that v minimises J\ in BV(fi) and seek height bounds. Note that 
a bound from above on v would correspond to a bound from below on u, while a strictly positive 
bound from below on v corresponds to a bound from above on u. 

At the end of this section we show an easier way to find one-sided estimates in the cases of j3 
being either positive (for which we show v bounded from above) or negative (for which we show v 
bounded from below). 



3.1. Height bounds from above on v. To estimate a minimiser v from above, we follow an ap- 
proach similar to [10, 15] leading to a Stampacchia iteration [18]. We use the following variant of 
the original Stampacchia lemma: 

Lemma 3.1. Suppose B(t), non-negative and non-increasing in t, satisfies 
(3.1) (h - k)B{h) < Ck [B(k)} 1 , 

for all h, k such that < ko < k < h, for some constants C, kg > and 7 > 1. Then B(K) = 
for some sufficiently large K dependent on C, 7, ko and B(ko). 

We will use the following in the proof of the above lemma: 



Lemma 3.2. For all a > 1 and d > —a, the sequence 



s m = 1 + 



d 



a" 



1 + 



d 



m—l 



a 



1 + 



a 



converges to a non-zero limit. 

Proof. We examine the sequence {In s n }, writing each term as the partial sum YljXi m (l + 
and using the ratio test for the convergence of series: 



lim 



In 1 + 



d 



lim 

j— >oo 

lim 

j~*QO 
1 

a 



9 , (, d 



9 1 ( , d 
— In H r 

dj V aJ 



-\ -1 



din a 

ai +1 + d 



din a 

a? + d 



-\ -1 



As this series converges to some limit L, {s m } converges to e L > 0. 



□ 



Proof of Lemma 3.1. We begin by defining the sequence k m := k s m , where s m is as in the pre- 
ceding lemma, with a = 2 and d = C [^(fco)] 7 " 1 > (we assume here that B(k ) / 0, 
otherwise the lemma is trivially true). Note that as k m+ i — k m = kod2~( m+1 ^ > 0, {k m } is positive 
and increasing, and, by the above result, converges to some limit K. 

We now prove that B(k m ) < B(k )2^ m for p, = (1 - 7)" 1 < 0, by induction. 

The base step, for m = 1, is as follows: by assumption (3.1), 



(fci - k )B(h) < Ck B(k y, 
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and so using k\ — ko = kod/2 we find that 

B(h) < C^B{k y- l B(k^ 

= 2 l -^l^- l ^B{k Q )2^ 
= B(k )2». 

Now we make the inductive assumption that B{k m ) < B(ko)2 fim . We use this and condition 
(3.1) to estimate 

orrt+l 

fl(Wi) < c—B{k m y 

om+l 

< [5(fco)2 Mm ] 7 

= 2 m+1 -~<l^-^B{k Q ) l -"i [B(k )2» m ]' y 

< B(k )2^ m+1 \ 

Finally, the monotonicity of B implies that B{K) < linim^oo B(k m ) < lim™-^ B(ko)2^ m = 



0. 



□ 



Theorem 3.3. Let v minimise J\ in BV + (f2), where J\ is given by (2.2). Assume in addition <9$7 
to be Lipschitz and to satisfy an interior sphere condition. Then v is bounded above. 



Proof. We set A(k) = {x G ft : v(x) > k} for k > ko, ko to be chosen later, the goal being to 
show that the non-increasing |^4(/c)| vanishes for some large k. We also write w := min(ti, k). As 
v minimises J\, we have J\(v) < J\(w) for all eligible w, which after rearranging gives 



(3.2) > 



\ [ ^/v 2 + \Dv\ 2 - [ kdx\+ [ [ In hdhdx+ [ p( 
[JA(k) JA{k) J JA(k)Jk JT, 



v — w) ds. 



Here, we make use of the fact that w G BV(fi), that Dw = Dv inQ \ A(k) and Dw = in A(k). 
We estimate the boundary term using (2.1): 



/ P(v-w) 

JT. 



ds 



< (1 - a) 



/ \D(v — w)\ + cr \v — w\ dx 
J Si JSi 



One can easily show that J Q \Dv\ < J Q ^/v 2 + \Dv\ 2 , so that (3.2) gives 

(3.3) k\A(k)\>a f \D(v -k)\ + (In k - (1 - a)c R ) f \v-k\dx. 

JA(k) JA(k) 
For BV functions on C 0,1 domains $7, we have the Sobolev inequality 



7 

JSi 



\f\ dx 



< C2 



\Df\+ / \f\dx 
Si JSi 



To see this, we first note that the inequality holds true for W l,1 (Vl) functions, since the space 
W 1 ' 1 ^) is continuously embedded in L n /(" -1 )(n) for n > 1 (see Theorem 7.26 in [11]). The 
results extends to / € BV(fi), after approximating / by smooth functions as in Theorem 1.17 of 
[13], and then following the steps of the proof of Theorem 1.28 of [13]. 
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We rearrange this inequality as 

n-l 

a [ \Df\>— f \f\^dx -a f \f\dx. 
Jn c n Un J Jn 

Using the above estimate with f = v — w, (3.3) becomes 



n-l 

CL f n n I" 

k\A(k)\ > — / \v — w\ dx + (In k — (1 — o)cr — a) I \v — w\dx. 
c n [Jn J JA{k) 

By choosing ko (the lower bound on k) large enough, we can ensure that (In k — (1 — a) cr — a) is 
positive, and drop this term altogether. The Holder inequality [13, Theorem 1.28] gives 

and so for all h, k with h > k > ko, we have 

— (h - k)\A(h)\ < — [ {v-k)dx< k\A(k)\ 1+ n. 
c n cn J A ( k ) 

Now we can apply Lemma 3.1 and conclude that for sufficiently large K, \A(K)\ =0. □ 

3.2. Height bounds on v from below. We start by remarking that bounds from below on v — 
or from above on u, where the surface is given by graph u — are not essential for the existence 
proof. We will see in the next section that the energy is bounded from below irrespective of such an 
estimate, and that subsequent results, leading to existence, also hold. However, they are important 
for the correspondance between u and v, and for the regularity results. 

Theorem 3.4. Suppose that v minimises J\ in BV + (f2). Then there exists a bound from below on 
v, < c < v, ^-almost everywhere. 

The proof is similar to that of Theorem 3.3, but we use a slightly stronger Stampacchia-type 
result. 

Lemma 3.5. Suppose B(t), non-negative and non-increasing in t, satisfies 

(3.4) (h-k)B(h) <Ch[B(k)}~ / , 

for all h, k such that < ko < k < h, for some constants C, ko > and 7 > 1. If 

(3.5) CB(ko)^ 1 < 1, 
then there exists a K < 00 such that B(K) = 0. 

Proof. From (3.5), we may choose d > and a > 1 such that CB(ko)' y ~ 1 a~ < d < a. 

Next, define the sequence k m := k /s m = fc m _i(l — ^-)~\ where s m is as in Lemma 3.2 (but 
note the change of sign on d); here {k m } is positive, increasing, and by Lemma 3.2, convergent to 
some K. 

— m 

We now prove that B{k m ) < B{ko)ct"<- 1 by induction. 

The base step, for m = 1, is as follows: by the definition of the sequence, we have k\ = 
fco(l — d/a)^ 1 , and so our assumption (3.4) gives 

(fci - ko)B{h) < ChB(koy, 
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which leads to 



Now we make the inductive assumption that B(k m ) < B{ko)a"<- 1 , and show that this then holds 
for krn+i- we use (3.4) to estimate 



b(Wi) < c - km+ \ B(k m y 



krn+1 k m 

a m+l 



iv -, —mf | m+1 

< C^— B(k ) 1 ~ 1 a— + — 

-(m + 1) 

< B(k )a 7-1 . 



-(m + 1) 

B(k )a 7-1 



Finally, the monotonicity of B implies that B (K) < linim^oo B(k m ) < linim^oo B(ko)a-i- 1 = 
0. □ 



We will need to show that the measure of the set where v is small is small enough to satisfy (3.5). 
This is the only place in this paper where it is not immediately clear that rescaling the constants 7$ 
to 1 does not result in a loss of generality. Consequently, we include the arbitrary weightings in J\ 
in the following step. 



Lemma 3.6. Let v minimise J\ = jiE$ + 72 W + 73-Es in BV + (f2), and set B(k) := {x <E 0, : 
v(x) < l/k}. Then for all rj > we can find a k such that \B(k)\ < rj. 



Proof. Define the comparison function w := maxju, M G BV + (f2) for any k > ko, ko to be 
chosen later. Note that < w — v < i. Since v minimises J\, we have J\(y) < J\(w). We use 
J \fu 2 + \Du\ 2 > J \u\ dx for u G BV + (0) to estimate 

o> J1W-J1W 

= 71 I yj v 2 + \Dv\ 2 — 71 I — dx — 73 J P(w — v)ds 

JB(k) JB(k) k JT, 



f , 1 , 1 1 , 

+ 72 / v Ixiv — v — — in - — I — — dx 

JB(k) 



> 



I -171 - 

JB{2k) 



(fc) 

72 1 



J(k) 

A 1111 A, 



A: 



1 



1> — 



B(k)\B(2k) 

- 73 (l-a)|S| 



k 

1 7i - 72 



+ "2 I r In r ^ In j 



dx 

+ 72 I v In v — 



A: 
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— now choose ko large enough so that x In x is decreasing for < x < ^ — 

>/ ~|7i -72|t +72 ( - 7: ln 7: J dx + I -|7i -72 1 ^rdx 

JB(2k) k \2k 2k k k) J B {k)\B{2k) 2fc 

- 73 (l-o)|E|i 

> -171 " 72|^| - 73d - fl )|E|i + 72 j ^ (-^) (In* + 1) dx 
for some e G (^,|) 

> -|7i - 12\\M ~ 73(1 - a)|S|i + 7215(2^)1 (-^) (in ±- + 1 
Rearranging, and choosing fco large enough that In ^ < —1, we find that 

|*(»)| < ^-nlPH-nd -'iH < , 

1 V ;I - 7 2 (ln2fc-l) 
for sufficiently large k > k . □ 

Proof of Theorem 3.4. Let B{k) be defined as above. Set w := max(ti, l/k), for some k > ko. 
Again, as v minimizes J\, then J\{v) < J\{w). Proceeding exactly as in the proof of Theorem 
3.3, we obtain 



> a [ \D(v -w)\- j\B(k)\ - [(1 - a)ci + ln(l/ife)] f 



\v — w\ dx 

n 



> — \\r - ir\\ , rii - ^\ B {k)\ - [(1 - a)ci + ln(l/fc) + a] ^ |u - io| dx, 



a 

~"^(^) fc 1 

and if we choose ko large, so that In A; > ci(l — a) + a, then the final term above is positive. 

n 

We drop it and apply the Holder inequality to the L«-i term, leaving us with (l/k)\B(k)\ > 
(a/cn)\\v — w\\ L i^ B ^\B(k)\~ 1 / n , and so for each h > k > ko we have 

(h — k)\B(h)\ < Ch\B(k)\ 1+ ^. 

Lemma 3.6 implies that we can find ko large enough that CB(ko)« < 1. We can then apply 
the Stampacchia-type Lemma 3.5 to conclude that \B(K)\ = for large K, and so v > almost 
everywhere. □ 



3.3. Height estimates in the cases (3 < and f3 > 0. Height estimates are easier to obtain in case 
(3 is either non-positive or non-negative. 

We begin by observing a height bound for v in the surface energy term. This closely follows 
Lemma 5 of [2], and may be proved in the same way. 

Lemma 3.7. Let v G BV + (f2) and suppose that A(k) = {i £ : v(x) > k} has positive 
measure. Then w = min(t>, k) £ BV + (f2) and for almost all k, 

[ ^w 2 + \Dw\ 2 < ( V« 2 + \Dv\ 2 . 
Jn Jn 
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Theorem 3.8. Suppose that v G BV + (f2) minimises J\, and that (3 > 0. Then v is bounded from 
above. 

Proof. Set w = min(v, k). Suppose that A(k) is of positive measure for some k > 1. We may 
choose A; so that Lemma 3.7 gives us 



/ ^w 2 + |L>w| 2 - / + |L>?;| 2 < 0. 
Jn Jn 



We note that 



In hdhdx — / In h dh dx = / / In hdhdx 
inJi JnJi JnJv 



-I I 

J nn{x:v(x)>k} Jv 

< 0. 



In /i d/t dx 



Finally, 

- / f3[l - w] ds + / p[l -v]ds= / /3[u> - v] ds < 
Js Je is 

if /3 > 0. Together, these inequalities give J7i — J7i (u) < 0, contradicting that v was a minimum. 
It follows that \A(k)\ cannot be positive, and so v < k. □ 

Lemma 3.9. Suppose that v minimizes J\, and f3 < 0. Then v > e~ l . 



Proof. Set w = max(t>, e), and write B(e) = {i £ !1 : f (x) < e}. Then 

r-max(t),e) 

V^ + I^l 2 " / / 

He) JB(s) 

Ids 



Ji(v)-Ji(w)= / yjv 2 + \Dv\ 2 - / V^ 2 + |£H 2 - / / Inhdhdx 

JB{e) JB{e) Jn Jv 

-J/3[w-v] 

= / \/v 2 + \Dv\ 2 — / edx— / Inhdhdx— (3[max(v , e) — v] ds 

Jb(e) Jb(e) JB(e) Jv Jt, 

> / ^Jv 2 + \Dv\ 2 - \ {e-v + v)dx + (-lne) / (e - v) dx 

•Ab( £ ) iB( £ ) ./£(e) 

= / \J v 2 + |Dt>| 2 — / vdx + (—lne — l) (e — v) dx 

Jb(e) Jb{e) Jb{e) 

> 

for all e < e _1 , if |e — v\ / 0. However, this would contradict our assumption that v is 
minimal for J\, so we conclude that \B[e)\ = for small enough e. □ 



4. Existence of a minimiser 

Lemma 4.1 (Lower bounds for the energy). If v G BV + (f2), then J\{v) > C(n, R, a, where 
C is not necessarily positive. 
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Proof. As before, we can incorporate the wetting energy into the surface tension term using (2.1), 
so that 

Ji(v)> f Jv 2 + \Dv\ 2 + [ [ Inhdhdx- (l-a)cji f |1 - v\ dx - (1 - a) f \Dv\ 
Jn Jn Ji Jn Jn 

> a / \Dv\ + / f(v)dx, 
Jn Jn 

where f(v) := v(\nv — 1) + 1 — cr(1 — a)|l — v\ is bounded below by a constant dependent on 
cr and a. The result follows. □ 

We define a minimising sequence for J\ as a sequence Vj € BV + (0) with 

lim J\{vj) = inf Ji{w) := m. 

Lemma 4.2. A minimising sequence for J\ is uniformly bounded in the BV-norm. 

Proof. We can assume that J\ (vj ) < m + 1 for j large enough. As in the previous lemma, where 
we defined /, we then have 

m + 1 > Ji{vj) > a \Dvj\ + / f(vj) dx, 
Jn Jn 

so the uniform bound follows from the lower bound on /: 

j a \D,,\ <I( m+ l -|fl| m K f + /(/,)). 

Also, since there exist positive constants ai,a2 such that fit) > ait — a-i, we have the uniform L 1 
bound 



'jhnn) < — { [ if(vj) +a 2 ] dx) < — (m + 1 + a 2 |^|). 



□ 



Lemma 4.3 (Lower semicontinuity of J{). A sequence Vk G BV + (f}) with Vk —>■ v in L 1 (J1) 
satisfies 

Jiiv) < lim inf Ji(v fc ). 

fe^oo 

Proo/ We show the surface energy term is lower semicontinuous. For any admissible g, we have 

/ v (g n+1 + div n g) dx = lim / v k (g n+1 + div n g) dx 
Jn k -^°°Jn 

= lim inf / v k (g n+1 + div n g) dx 

< lim inf / Jv\ + \Dv k \ 2 . 
k ^°° Jn v 

Lower semicontinuity follows by taking the supremum over all admissible g. 

Continuity of the remaining terms of J\ follows as in [10, Appendix II]. □ 

Combining all of the above results we have: 

Theorem 4.4 (Existence of a minimiser). There exists a function v € BV + (f2), such that 

Mv) = inf J x {w). 

w&BV+{n) 
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Proof. Let {vj} be the minimising sequence of Lemma 4.2, with \\vj ||bv(o) < C. By the standard 
compactness theorem (for example [13, Theorem 1.19]) there exists a subsequence Vf — > v in 
L 1 (J7). 

Since the BV-norm is lower semicontinuous, f is also in BV + (f2), andj^u) > inf 1)gBV +(Q- ) i7i(u;) 
m. Lower semicontinuity of J7i, as in Lemma 4.3, gives J\(v) < liminf J\(yji) = m, completing 
the proof. □ 

5. Regularity 

In this section we show that a minimiser v G BV + (f2) of J\ is locally Lipschitz in 17 following 
a procedure similar to [10]. 

In a subsequent paper, we discuss boundary regularity. If one has boundary regularity, the meth- 
ods of [1] can be used to derive higher regularity in smooth domains. 

Theorem 5.1. Let v be a minimiser of J\ in BV + (17). Then v is locally Lipschitz in fi. 

Proof. We mollify v over the whole of Q. The mollification v £ is in C°°(Q), and shares the height 
bounds derived for v in Section 3 (that is, bounded above and bounded from below away from zero). 
Furthermore, since v G BV(fi), 



v £ — > v in and / \Dv £ \ — > / \Dv\ 

Jn Jn 



Corollary 2.3 for the surface energy and standard convergence results for the remaining energy 
terms then imply that 

(5.1) Mv £ )^Mv). 

Let B C VL be any ball of sufficiently small radius p, and consider the following two related 
Dirichlet problems: 

y/w e 2 + \Dw e \ 2 ^/w £ 2 + \Dw £ \ 2 

w £ = v £ on dB; 



div . f Wg = = . n w - + In w £ in B, 
(5.2) ....... ,■ 



and 
(5-3) 



div . ^ o = a _J = + u £ \ in B, 



ti e = — cringe) on 

The second expression is in fact a family of problems, indexed by a £ [0,1]. This family is of 
mean curvature type. Note that for smooth w £ and u £ , (5.2) is equivalent to (5.3) for a = 1, with 
the correspondence w £ = e~ Ue . 

Our next step is to solve (5.3) for a = 1 using the continuity method. We apply [11, Theorem 
13.8]. A prerequisite for this is to show that a smooth solution u u of (5.3), for any a G [0, 1], has 
height and gradient bounds independent of a. 

The height bound may be found in [17]; however, the geometric nature of our problem admits a 
shorter proof which we present as the following lemma. 
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Lemma 5.2. Let u a be a smooth solution to (5.3) corresponding to a a £ [0, 1]. Then 
(5.4) sup|u CT | < Mi, 

B 

where M\ depends only on sup^ |lnw e |. 

Proof. We suppose that u a achieves a positive interior maximum, u a (x) = M at some point x £ B. 
If M > 1, then the mean curvature H(u a ) = div , Du " = at x must be strictly positive. But 
a point of positive mean curvature cannot correspond to an interior maximum, contradicting the 
assumption M > 1. We conclude that 

u a < min{l,sup |oTnv e |} < min{l,sup | lnw £ |}. 

dB ' dB 

A similar argument shows that u a has no negative internal minimum, so u a > — sup^s | In v £ \ . □ 



Continuing the proof of Theorem 5.1: We find that the gradient bound 

(5.5) sup\Du a \<M 2 

B 

is an application of standard results. Firstly, an interior gradient bound can either be derived by 
applying a maximum principle to the elliptic equation satisfied by the gradient; or by using [16, 
Theorem 4], which gives 

(5.6) sup \Du a \ < M 3 

B> 

where B' CC B and M3 is dependent on dist(f?', dB), n and sup \ u a \. 
Secondly, a boundary gradient estimate 

sup|Du CT | < M 4 

dB 

results from [11, Corollary 14.5] with the structure condition (14.33). Here M4 is dependent on 
I h\v £ \ C 2( dB y n, sup \u a \, and p. Together these two gradient estimates give us (5.5). 

The conditions for the continuity method being satisfied, the problem (5.3), with a = 1, has a 
C 2 ' a (B) solution which we call u £ . It has height and gradient bounds (5.4) and (5.5). It is also 
unique: the proof is similar to that of Theorem 2.2 in [1], adjusted to Dirichlet boundary data. 

We set w £ = e- u ^. This is a C 2 > a (B) solution of (5.2) with height bound e~ Ml < w £ < e Ml 
and gradient bound \Dw £ \ < M2e Ml . 

Note that (5.2) is the Euler-Lagrange equation for the energy 

J 2 (w):= J ^w 2 + \Dw\ 2 + J Inhdhdx, 

and so w £ is a critical point of J2 in the class of H 1,2 (B) functions with boundary data v £ . Fur- 
thermore, as the integrand of J2 is convex in (w, Dw), w £ is also a minimiser in this class (see, for 
example, the remark in Section 8.2.3 of [5]) and hence in the smaller set C 2 ' a (B). 
In particular, if we compare w £ to v £ , we have 



(5.7) J ^w 2 + \Dw £ \ 2 + j £ Inhdhdx < J y ' v 2 + \Dv £ \ 2 + J Inhdhdx. 
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Now let v £ be denned by 

w £ in B 
v £ in $7 \ B. 

Using (5.7) for the region B where v £ may be different to v £ , we see that J\{v £ ) < Ji{v £ ). 

Now we will show that v £ converges to a BV(fJ) function which is locally Lipschitz. 

Uniform L 1 (fi \ B) bounds are given by the height bounds for v in Section 3. Uniform L 1 (-B) 
bounds are given by sup B \ w E \ < e Ml where Mi is the constant in (5.4); Mi also depends on the 
height bounds for v. 

As a consequence of (5.1), we may assume that Ji{v £ ) < Ji{v) + 1. Then J\(v £ ) < Ji(v) + 1, 
and so 

[ Jvl + \Dv £ \ 2 < Ji{v) + 1 - [ [ \nhdhdx+ [ p[l-v £ ]ds 
Jn Jn Ji Jt, 

<Ji(v) + l + \n\ sup (hlnh - h+ 1) + |S|sup|/3|(l +sup|w|) 

inf v<h<supv S f2 

which is bounded above, independently of e. Uniform BV bounds follow as in Lemma 2.2. There- 
fore a subsequence of v £ converges to v q € BV(S7), and vq is Lipschitz in B' with bounds given by 
(5.6). 

Lower semicontinuity of the functional now gives 

Ji(«o) < liminf J x (y £ ) < liminf J ± (v £ ) = Ji{v) 

but as v was assumed to minimise Ji these must all be equal. We conclude that there exists a 
minimiser of Ji that is locally Lipschitz on interior sets. □ 

Reconsidering the problem of a capillary surface S = graph u that minimises the original energy 
functional J given in the introduction, we conclude that the found v corresponds to a minimiser in 
the class {w : e~ w G BV + (f2)}. This solution is given by u = — \nv, and is locally Lipschitz on 
interior sets. 
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